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Abstract: The energy evolution of medium-modified average multiplicities and multiplicity fluctuations 
in quark and gluon jets produced in heavy-ion collisions is investigated from a toy QCD-inspired model. 
In this model, we use modified splitting functions accounting for medium-enhanced radiation of gluons 
by a fast parton which propagates through the quark gluon plasma. The leading contribution of the 
standard production of soft hadrons is found to be enhanced by the factor \/]V7 while next-to-leading 
order (NLO) corrections are suppressed by l/^/Wg, where the nuclear parameter A'^s > 1 accounts 
for the induced-soft gluons in the hot medium. The role of next-to-next-to-leading order corrections 
(NNLO) is studied and the large amount of medium-induced soft gluons is found to drastically affect 
the convergence of the perturbative series. Our results for such global observables are cross-checked and 
compared with their limits in the vacuum and a new method for solving the second multiplicity correlator 
evolution equations is proposed. 
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1 Introduction 



The properties of quark and gluon jets are strongly established and carefully studied through well known 
QCD evolution equations [1] in e^e~, ep and hadron colUsion experiments. In the case of high en- 
ergy nucleus-nucleus collisions, hard jets propagate in a medium with different properties from those 
of the vacuum. Recent experiments at the Relativistic Heavy Ion ColUder (RHIC) have established a 
phenomenon of strong high-transverse momentum hadron suppression [2, 3], which supports the picture 
that hard partons going through dense matter suffer a significant energy loss prior to hadronization in the 
vacuum (for recent review see [4]). 

Since little is known so far on jet evolution in QCD media, predictions for multiparticle production in 
such reactions could be carried out by using a toy QCD-inspired model introduced by Borghini and 
Wiedemann in [5]; it allows for analytical computations and may capture some important features 
of a more complete QCD description. In this model, the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi 
(DGLAP) spUtting functions q ^ gq and g ^ gg [1] of the QCD evolution equations were distorted so 
that the role of soft emissions was enhanced by multiplying the infra-red singular terms by the medium 
factor Ng: 



where x is the fraction of the outgoing jet energy carried away by a single gluon. Thus, the lead- 
ing singular terms of the spUtting functions play a more important role; from the theoretical point of 
view it could be considered as a result of some effective Lagrangian, which would be responsible for 
processes in a dense nuclear environment. In the approach [6], the medium-modified splitting func- 
tions are directly related to the medium-induced gluon spectrum dl'^^/dxdE [8], where E is the initial 
energy of the emitting parton going through the medium. As compared to the Borghini-Wiedemann 
model, the medium modifications explicitly depend on the parton virtuality through the enhanced in- 
duced gluon spectrum. However, we use the simpler interpretation of the induced-medium modification 
of the Borghini-Wiedemann model [5], which was further discussed and used on the description of final 
state hadrons produced in heavy-ion collisions [7]. 

In this paper we are concerned with multiparticle production in quark and gluon jets A = q,g, pro- 
duced in nucleus-nucleus collisions at very high energy. We make predictions at NLO and NNLO for 
the average multiplicities Na, for the ratio r = Ng/Ng and finally for the second multiplicity correlators 
{Na{Na — 1))IN\, which defines the width of the multiphcity distribution. Such observables are of 
great importance and have been largely studied in the vacuum from both theoretical [10-13] and experi- 
mental [14—19] points of view. The problem of medium-modified multiparticle production has also been 
considered in [21,22] with fixed coupUng constant. 

The starting point of our analysis is the NLO or Modified-Leading-Logarithmic- Approximation (MLLA) 
evolution equations [1], which determine the jet properties at all energies together with the initial condi- 
tions at threshold at small x. Their solutions with medium-modified sphtting functions can be resummed 
in powers of ^J ag/Ng and the leading contribution can be represented as an exponential of the medium- 
modified anomalous dimension, which takes the A'^s-dependence: 




Na — exp 
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where 7med(«s) can be expressed as a power series of ^fajWl in the symbohc form: 

7med {as) ^ X ^ (^1 + + ^ + . . 

Within this logic, the leading double logarithmic approximation (DLA, 0{\/l%a^)), which resums both 
soft and coUinear gluons, and NLO (MLLA, 0{as)), which resums hard collinear partons and accounts 
for the running of the coupling constant Ug, are complete. The DLA takes into account, as expected, 
the medium modification by enhancing the soft multiparticle production by a factor oc the MLLA 
terms, which are A'^s -independent, takes into account the energy balance in the hard collinear parton split- 
ting region as in the absence of the nuclear modification. However, this result depends on the definition of 
the medium-modified splitting functions. The next terms, which are NNLO or next-to-MLLA (NMLLA, 
0(05^^/^/^)) are not complete but they include an important contribution, which takes into account 
energy conservation and provide an improved behavior near threshold. With medium modification, the 
NMLLA terms take A'^g-dependence, but this will be explained in the main core of the paper. This logic 
applies to each vertex of the cascade and the solution represents the fact that successive and independent 
partonic splittings inside the shower, which in this case concern both vacuum and medium-induced soft 
gluons, exponentiate with respect to the evolution-time variable Y {dY = dQ/Q), where ^ 1 is the 
angle between outgoing couples of partons. The choice of y ~ ln(0) follows from Angular Ordering 
(AO) in intra-jet cascades; furthermore, the tree Feynman diagrams describing the process are at the heart 
of the parton shower picture [1]. Thus, the solutions of the equations incorporate the Markov chains of 
sequential angular ordered decays and 7med determines, in this case, the rate of multiparticle production 
in the dense medium. 

At the end of the cascading process inside the medium, partons hadronize in the vacuum. In order to 
obtain the hadronic spectra, we advocate for the Local Parton Hadron Duality (LPHD) hypothesis [23] : 
global and differential partonic observables can be normalized to the corresponding charged hadronic 
observables via a certain constant /C that can be fitted to the data, i.e. N\ = JC x Na- 

The paper is organized as follows: 

• Section |2] presents a system of evolution equations with medium-modified splitting functions, 
which allows for the computation of the medium-modified average multiplicity and the medium- 
modified gluon to quark average multiplicity ratio at NLO and NNLO. We give predictions for the 
values Ns = 1.6 and Ns = 1.8, which may be realistic for RHIC and LHC phenomenology [5]. 
Moreover, we compare our results with previous predictions in the vacuum; 

• in Section [3] we study the medium-modified second multiplicity correlator at NLO and NNLO. 
Accordingly, we give predictions for the same values of Ng and compare such predictions with the 
equivalent for the vacuum limit A'^s = 1; 

• in Section|4]we present our conclusions. 
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2 Evolution of the average multiplicity and gluon to quark average mul- 
tiplicity ratio with energy loss 



At MLLA the evolution of quantities with jet energy E and jet opening angle 6 is given by an evolution 
equation for the azimuthally averaged generating functional in the jet [1]. The evolution involves Og, the 
running coupUng constant of QCD: 

as = asiQ) = ^^^TTV' ^0 = ( - ^Tr] , (1) 



4iV,/3oln(S) 



4iVc V 3 3" 



where Q = EQ is the maximum transverse momentum of the jet, A = Aqcd is the intrinsic scale of 
QCD, Po is the first term in the perturbative expansion of the /?— function, Nc is the number of colors and 
Tr = nf/2, where nj is the number of quark flavors. In the leading DLA, ag is linked to the anomalous 
dimension 70 of twist-2 operators by the formula: 

7o = 7o(Q) = 2A^e^ = ^-^^, ^ = 1-Q^> A = ln-, (2) 

where Qo is the coUinear cut-off parameter kx = EQ > Qq. The results depend on energy and angle 
only through the variable Y. For the sake of simplicity we also set y' = y + A in the following. The 
average multiplicity is obtained by integrating the one-particle single differential inclusive cross section 
over the energy fraction x = e/E 



For the medium-modified evolution of the average multipUcity in quark and gluon jets one obtains as a 
consequence of AO at MLLA, the coupled system of two evolution equations 

^^9iY) dxjiix) [^liN,, x) {Ng{Y + lnx)+ Ng{Y + ln(l - x)) - Ng{Y)) 

+nf^l{Ns,x) {Ng{Y + lnx) + Ng{Y + ln{l - x)) - Ng{Y))] , (3) 



_d_ 
dY 



NgiY) = f dx^l{x)[mN,,x) {Ng{Y + \nx) + Nq{Y + ln(l - x)) - Nq{Y))\ (4) 



with medium-modified splitting functions as suggested in [5] in the Borghini- Wiedemann model 

$^(iV„x) = ^ - (1 - x)[2 - x{l - x)l ^l{Ns,x) = + (1 - x)% (5) 

which accounts for energy loss in the medium by enhancing the singular terms like $ Ng/x as x ^ 0. 
The g ^ qq spUtting function as well as the regular parts of g ^ gg and q ^ gq splitting functions are 
hard and provide colUnear corrections, that is why these terms do not take Ng dependence. 
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2.1 MLLA evolution of the average multiplicity in the medium 



For y » Inx ~ ln(l - x), N{Y + Inx) {N{Y + ln(l - x))) can be replaced by N{Y) in the hard 
partonic splitting region x ~ 1 (1 — x ~ 1) (non-singular or regular parts of the splitting functions), 
while the dependence on In x is kept on the singular one ^{x) Ng/x as it is performed in the vacuum. 
Furthermore, the integration over x can be replaced by the integration over Y(x) = In (^^^^^ 

7o'(^)- = / li{Y{x))dY{x). (7) 



After applying the differential operator ^ to both members of the system ( I3I4I ) above, one is left with 
the approximate system of coupled equations, 

NgiY) = ^nN,-m—)Ng(Y), (8) 



dY^ ' "dY 



N,{Y) = ^7o' {ns - ai^) N,{Y), (9) 



dY"^' 

with the initial conditions at threshold A^^(O) = 1 and A'^^'^(O) = and the hard constants 



ai 



3 3 I Nr 



a, = -. 



The quantum corrections oc ai , ai in ( I8I9I ) arise from the integration over the regular part of the splitting 
functions, they are ©(-y/a^) suppressed and partially account for energy conservation as happens in 
the absence of the dense medium. Since only the DLA terms are medium-enhanced in (I8I9I) . the hard 
constants are A^^ -independent. 

These equations can be solved by applying the inverse Mellin transform: 

to the self-contained gluonic equation where the contour C lies to the right of all singularities of 
Ng{i^) in the complex plane. The running of the coupling constant agiY), Eq. Q, is taken into account 
through the identity 

Consequently, one is left with the following differential equation in Mellin space 

NJuj) = ^+ — -2 -. (11) 

Solving (fTTl) and using (fTOl) . one obtains 

The exact solution of (fT2l ) together with the initial conditions leads to a linear combination of two kinds 
of Bessel functions which resums the perturbative series at all powers of [1]. However, in this 
paper we are concerned with the asymptotic solution of the equation as y » 1 {EQ ^ Qo), that is 
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the high energy Umit. Therefore, the MeUin transform (fT2b can be estimated by the steepest descent 
method. Indeed, the large parameter is Y' and the function in the exponent presents a saddle point at 
ujQ = ^Ng/iPoY'), such that the asymptotic solution reads 



n!'(Y) ~ /C X y 1^0 exp . 



/3o 



-Y' 



(13) 



where 

We also introduced, as stressed in the introduction, the LPHD normalization constant /C [23], which 
accounts for hadronization effects outside the medium. The constant ai is A^^ -independent because 
it resums vacuum corrections. Therefore, the production of soft gluons in a dense medium becomes 
2( 



exp 



/n:-i)^/yWo 



higher than the standard production of soft gluons in the vacuum [1] and the 
factor y/W^ in ([T3] ) underlines the presence of the nuclear medium; this results has first been reported 
in [24]. From ( [T3l ) one obtains the first and second logarithmic derivatives of Ng-. 



1 dNg 

WgllY 



1 dNg 
'Ng'dY 



1 d^Ng 
Wg~dY^ 



1 d'^Nq 



(14) 



The expression on the left hand side of ([T4l) is nothing but the MLLA rate of multiparticle production 
with respect to the evolution-time variable Y ~ ln(0) in the medium, which we define as the medium- 
modified MLLA anomalous dimension: 



7 d(y) = -^ 



1-ai 



70 



+ 



(15) 



where Ng only affects, as expected, the leading double logarithmic term. From ([T3] ) and ([T5] ). one 
recovers the ansatz 

Ng{Y) ~ exp f / i^,A{Y)dY\ , (16) 



which we further improve in the next paragraph by adding higher order corrections. Finally, using ([13 
and (I8I9I ). one obtains the solution for Nq-. 



N,{Y) 



Cf 

Nr. 



1 + (ai - di] 



70 



NgiY) + 



N, 



(17) 



Therefore, we can introduce the medium-modified MLLA gluon to quark average multiplicity ratio r 

Ng/Ng = Nj^/Nj^ in the form 



r = ro 



WK \n. 



N^ 
Cf' 



ri = ai - ai, 



(18) 



where the suppression factor restricts the production of hard coUinear partons as A'^s > 1- We 

notice that ([TSl l is identical to the expression with fixed coupling constant as{Y), where F = In (Q/Qo) 
and Q = EQ is the virtuality of the jet produced in the nucleus-nucleus reaction. This factor is found to 
suppress the hard correction O(7o) and therefore, r approaches its asymptotic DLA limit tq = Nc/Cp 
when the coherent radiation of soft gluons is enhanced. Setting A^s = 1 in (fTSl) . one recovers the 
appropriate limits in the vacuum [1,1 1,25]. The constants entering in (IT3l) and (fTSl) are the same as those 
obtained in the vacuum and their values are displayed in Table [U 
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Uf ri 0"! 



3 0.185 0.280 

4 0.191 0.297 

5 0.198 0.314 



Table 1 : Coefficients ri and cji . 



2.2 Medium-modified equations and solutions at Next-to-MLLA 

Previous MLLA results for the average multiplicities can be improved by further pushing the perturbative 
series in ( 1314 1 ). We can include NNLO or NMLLA corrections of order 0{as), which are known to better 
account for energy conservation in the vacuum [12, 17,20]. For this purpose, we replace N{Y + Inx) 
(iV(y+ln(l-x)) by the Taylor expansion N{Y) + -^N{Y) lnx+. . . {N{Y) + -^N{Y) ln(l-x) + . . .) 
and, as in section 12.11 we integrate over the non-singular parts of the splitting functions. We thus obtain 
the medium-modified NMLLA approximate system of two-coupled evolution equations 

^^.(^) = 7o' {ns - «i - f3li) + «2(iV.)^) N,{Y), (19) 
N,iY) = %l (n, - ai f a _ a^l] + ~a,{N,)4^\ N,{Y), (20) 



with the new A'^s -dependent constants 

(AT^ 67 7r2 UTrCf , 2TRCF {ai-~ai) . , . 7 CpfB vr^ 

^2 A, = A, ^= — , ao A, = — A, — 

^ ^ 36 6 18 Ac Ac 3 Ac Ac ^/Ws 8 Ac V8 6 

The dependence of 02 and 02 on A^ follows from the singular term in the integral A^ ^ ln(l — x) = 

2 

— A^^ which also enhances (see below) the induced soft gluon radiation at the NNLO level. The term 
oc 1/\/As in 02 (As) was obtained by replacing ([TT] ) in the single logarithmic piece oc Ng in while 
the term oc A^ in the equations d 191201 ) enhances the role of the leading DLA as in ( I8I9I) . As expected, 
one recovers the constants 02 (A^ = 1) = a2 and 02 (A^ = 1) = a2 obtained in the vacuum [20] when 
As = 1. The terms proportional to these constants are known to better account for energy conservation 
in the partonic shower in the vacuum. Also note that the contributions oc aiPo and ai/3o cannot be 
neglected when performing predictions with running coupling coupling constant. 

The system (I19I20I ) can be solved by inserting the ansatz ([T6l) in both sides of ( [191 ) with 



JmediY) = ^/Ns-fo 



/As ^^As ' \ 



(21) 



where (T2(As) is the unknown coefficient to be determined. The medium-modified NMLLA anomalous 
dimension (l2Ti ) has been inspired from the MLLA ( [T5] ) where, in both cases, we make appear the rescal- 
ing of the coupling constant Os ois/Ng in the series. After equating terms oc 7q in the left and right 
hand sides of ([T9l ) we obtain the following value 

fT2(As) = - i Qai/3o + \al + Asa2(As) + ^/J^) , (22) 



6 



such that after integrating (|2T]) according to ( fT6l ). the medium-modified NMLLA average multiplicity 
takes the simple form 

^4iV^ 2(72 (iVs) 



exp 



'-Y' + 



(23) 



The term oc c72(-^'^s) in (l23l) provides the NMLLA correction to A'^q. As Ns increases, cr2(A'^s) follows 
the leading behavior a2{Ns) ~ A'^^ (and the second factor in the exponent of (l23l) ~ Ns^"^) which 
enhances the production of induced soft gluons in the medium at the NNLO level. Setting A'^s = 1 in 



(1231 ) one recovers the value of this constant in the vacuum cj2(A'^s = 1) = —0.378 as given in [12]. 
Furthermore, setting /3o = in (|22] ). the appropriate limit with fixed coupling constant = const can 
be deduced. 

We proceed to determine the medium-modified NMLLA gluon to quark average multiplicity ratio by 
subtracting (|20l ) from ( [19] ). one has 



(24) 



Applying the operators and to (|20|) and dropping corrections contributing beyond NMLLA, we 



obtain respectively 



~ Cf dY^ + "'^0 dY^^ ' ^""^ dY^ Cp dY^ 
where we keep track of the nuclear factor A'^s and the running of as(y(x)) in such a way that 
be rewritten in the form 



,d^Ng 



N, d'^N, 



(25) 



can 



A^„ 



Up Up 



Nc (ai-ai) dNg Nc /ai(ai-ai) (a2(A^s)-a2(iVs) A d^Ng 



after reassembling terms oc and oc Using ([14)) together with the initial conditions at threshold 
yields for r = Ng/Nq with NNLO accuracy 



dY 

cfNa 



Cp 



■+- 



dY^ ' 



70^ 



A^: 



3/2 



(26) 



where the coefficient r2 is explicitly dependent on Ns through the formula 

r-2(iVs) = (ai - ri + (a2(Ar,) - a2{Ns))Ns. (27) 

The dependence of r2{Ns) on [3q in ([26] ) underlines the account of the running coupling constant, setting 
/3o = the appropriate fixed coupling solution can be deduced. For A'^s = 1, the appropriate limits in the 
vacuum are recovered [12], for example 

r2{Ns = 1) = (ai - (Ti) n + (a2(l) - 02(1)) ^ 0.426. 

We give the NMLLA coefficients a2{Ns) and r2{Ns) defined in ([23] ) and ([26] ) as a function of A'^s in 
Tabled 

As Ng increases, the C'(7o/\/i^) correction decreases, while the one 0{'^q/Ns) becomes sizable and 
decreases Uke ~ —Ng. That is why it might be wondered whether the convergence of the perturbative 
series could be reached at a certain level of accuracy. Since the series widely oscillate at low energy 
scales, large terms oc vr^ in a2{Ns) and d2{Ns) might spoil or drastically affect the trends obtained at 
MLLA. This kind of behavior has first been noticed in the Koba-Nielsen-Olsen (KNO) problem [28] in 
the vacuum. 
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3 


0.087 - 




- 0.548iV^ + 0.914iV2 


-0.337 


- 0.014\ 


/n; 


- 0.850iVs + 0.822iV2 


4 


0.087 - 




- 0.494iV^ + 0.914iV2 


-0.319 


-0.019\ 


/n; 


- 0.823iV3 + 0.822iV2 


5 


0.087 - 


0.049V^ 


- 0.441iV, + 0.914iV2 


-0.302 


- 0.024\ 




- 0.797iV, + 0.822iV2 



Table 2: Coefficients r2{Ns) and a2{Ns). 
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Figure 1 : MLLA ([T3l) and NMLLA (1231 ) average multiplicity as a function of Q = EQ in the vacuum 
{Ns = 1) and in the medium {Ns = 1.6 and N., = 1.8). 



2.3 NLO and NNLO results on and r 

9 

Setting iif = 3, we display in Fig. [T] our results for the medium-modified MLLA ( fT3] ) and NMLLA 
(1231 ) average multiplicity. We plot Ng in the range 10 < QiGdV) < 500, where Q = EQ is the total 
virtuality of the jet related to Y in Q- We compare our results in the medium for Ns = 1.6 and A'^s = 1.8 
(see [5]) with predictions in the vacuum (Ns = 1), we set Qq = Aqcd = 0.23 GeV in the limiting 
spectrum approximation [12], and JC = 0.2 is taken from [12]. The values A'^s = 1-6 and Ns = 1.8 
in the medium may be realistic for RHIC and LHC phenomenology [5, 26]; the jet energy subrange 
10 < (3(GeV) < 50 displayed in Fig.[T]has been recently considered by the STAR collaboration, which 
reported the first measurements of charged hadrons and particle-identified fragmentation functions from 
p-i-p collisions [27] at ^snn = 200 GeV. Finally, the jet energy range in the same figure, in particular for 
those values at Q > 50 GeV, will be reached at the LHC, i.e Q = 100 GeV is an accessible value in this 
experiment (see [5] and references therein). 

Notice that at NMLLA, the increase of N^ with Ng is more substantial than at MLLA. The former is 
driven by the leading contribution to a2{Ns): it increases like (T2{Ns) ~ Ng (see Tabled in the sub- 
leading piece of (1231 ). In both resummations schemes we find, as expected from our calculations, that the 
production of soft hadrons increases as Ns > 1, which implies that the available phase space becomes 
restricted for the production of harder coUinear hadrons. 
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Figure 2: MLLA (fTSl ) and NMLLA (l26l) gluon to quark average multiplicity ratio as a function of Q 
E@ in the vacuum (Ng = 1) and in the medium {Ns = 1.6 and Ns = 1.8). 



In Fig. [21 we display the medium-modified MLLA ratio r = Ng/Nq (fTSl) and the medium-modified 
NMLLA ratio as a function of Q = EQ. As expected from (ITSl) . as Ng increases, the 0(70) correction 
is suppressed by l/^/Wl, the ratio approaches the DLA asymptotic regime rg = Nc/Cp = 9/4. At 
NMLLA, the previous trend goes in the opposite direction: one has indeed —r2{Ns) ~ —N^ (see 
Table 121) in the 0{'^q/Ns) piece of (|26l ). which is negative and sizable as Ng increases and, therefore, 
spoils the behavior obtained at MLLA signaling difficulties with perturbative theory in the medium. 
Indeed, sizable oscillations have been noticed in the perturbative series [13] and it turns out that they 
are wider in the medium than in the vacuum. For example, the NMLLA correction to r is 10% 
for Ng = 1, and for Ns = 1.8 it is ~ 40%. It should be noticed that every logarithmic derivative 
of Nj^ provides a half power of Ns to successive terms in the series in the form [NsUsY^^'^ 
(n = 1 ^ MLLA, n = 2 ^ NMLLA. . . ), such that the perturbative approach should fail as higher order 
terms are incorporated. Former statements suggest that by incorporating NMLLA and next-to-NMLLA 
(NNMLLA, see paragraph |3.4b corrections on an equal footing, the MLLA behavior can be recovered. 
We conclude from this analysis that MLLA provides a more realistic physical picture of the softening of 
jets than NMLLA. Therefore, either the incorporation of NNMLLA terms or the exact numerical solution 
of the evolution equations [7,29], which exactly accounts for the running of as and the energy balance, 
is required. In [21], a numerical solution of the equations was provided with fixed coupling constant a^, 
and the results are shown to follow our MLLA expectations as Ns increases. 

Finally, in both MLLA and NMLLA, the gluon jets are still more active than the quark jets in producing 
secondary particles and the shape of the curves are roughly identical; however, these characteristics prove 
not to be very sensitive to A^^- 



3 Medium-modified evolution for the second multiplicity correlator 

The second multiplicity correlator was first considered in [11] at MLLA and later in [13] at NMLLA. 
It is defined in the form N\ = {Na{Na — 1)) in gluon (A = g) and quark {A = q) jets. The 
normalized second multiplicity correlator defines the width of the multiplicity distribution and is related 
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to its dispersion D\ = {Na{Na — 1)) — N'\ by the formula 



n2 



{A2 - l)Nl + Na. 



(28) 



The second multiplicity correlators normalized to their own squared average multiplicity are 



Go 



{Ng{Ng-l)) 



iV2 
9 



Q2 



{N,iN, - 1)) 



Ar2 
g 



(29) 



inside a gluon and a quark jet respectively. These observables are obtained by integrating the double 
differential inclusive cross section over the energy fractions xi = ei/E and X2 = 62 /E 



{NAiNA - 1)) 



dxidx2 



1 (f 



a 



a dxidx2 J A 



The correlators are /C-independent and provide a pure test of multiparticle dynamics. However, important 
disagreements with e^e^ data [16] indicates that non-trivial hadronization effects may play a significant 
role. The treatment of this observable with full account of perturbative and non-perturbative ingredients 
is not available yet. Therefore, we study the variation of this observable as A^^ > 1 with respect to the 
limit in the vacuum A^^ = 1 and do not compare our results with e^e~ data [16]. The medium-modified 
system of two-coupled evolution equations for this observable follows from the MLLA master equation 
for the azimuthally averaged generating functional and can be written in the convenient form 



d 

dY 



(Ar(2) 



Nf\Y + \nx) + [Nf\Y + \n{l-x))-N^^\Y)^ 



+nfl dx-fi<^l 



+ {NgiY + lnx)- NgiY)) [NgiY + ln(l - x)) - Ng{Y) 



n^'\y)-n^{y) 



Y + \nx)-N^{Y + lnx 



+ {2Nq{Y + Inx) -Ng{Y)) (2Nq{Y + ln(l - x)) - Ng(Y) 



(30) 



dY 



(Ar(2) 



dxj^M [ivf )(y + inx)+('ivf )(y + in(i - x)) - ivf ^(y) 



+2[Ng{Y + Inx) - N,{Y) Ar^(y + ln(l - x)) - N,{Y) 



(31) 



which proves to be more suitable for obtaining analytical solutions in the following. We use a new 
method to compute solutions at MLLA and NMLLA by replacing N^^ = 7^2 -^1 on both sides of the 
expanded equations atx~ 1 — x~ 1. This observable is less inclusive than the average multiplicity, 
it can indeed be derived from the two-particle four-momentum correlation [30] in the shower. The 
medium-modified formulae (IT6l ). ([TSl l and (l26l ) will be used in this analysis. 



3.1 MLLA approximation 

For y > In X ~ ln(l - x) in the system above (I30I31I ). N{Y + In x) {N{Y + ln(l - x))) and iV^^) (y + 
Inx) (A^(^)(y + ln(l — x))) can be replaced by N{Y) and N^'^\Y) respectively in the hard partonic 
splitting region x ~ 1 — x ~ 1, while the dependence on Inx is kept on the singular one (x 0). 
The medium-modified MLLA approximate system of two-coupled evolution equations for the second 
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multiplicity correlator reads 



dY ) 



dY^ 

where the new hard constant is: 



iV{2) _ ^r2) = ^ 2 ~ A\ ^(2) 



ai— I ^v-' + (ai - ^i)7o^^g 



(32) 
(33) 



bi 



4iV. 



11 



—Nc - 4— 1 

3 A^. 



The constant Ng only affects the leading double logarithmic term of the equations. The terms propor- 
tional to oi, (oi — 6i) and di are hard vacuum corrections which partially account for energy conser- 
vation, indeed « a/^7o ^^'^ the relative correction to DLA is 0{^y ag/Ng). As before, the hard 
constants are, as expected, A'^s -independent. Moreover, these sub-leading contributions have the same 
form as those describing the two-particle correlation [30]. 



3.1.1 Medium-modified G2 at MLLA and expansion in 0(70 / V Ng) 

Setting A'^^^^ = G2Ng in (l32l) . the system can be solved iteratively by considering terms up to 0{a^^'^) 
in the left and right hand sides of (l32l ). Thus, the l.h.s. of (l32l ) writes in the form 



d^ 



ATf - AT^) 



d^G2^,2 dG2 d . 

+ 2 — A^^ 

dY^ 3^ dY dY <^ 



N + {G2 - 1) 



^2 



(34) 



dY^\ 3 9 J 9 ' dY dY ' " 'dy2 a- 

Hereafter, in all sub-leading terms, we can replace G2 by a constant G2 = G^^^ = const, while the 
terms involving Ng should be computed by using ([T3] ). thus 

d^ 



dY"^ 

while the rh.s. reads 



{Nf^-Nt)=llNs{G2 



1 



(4ai -/3o) 



70 



Nsll - "1^0;^) ^2 + (ai - ^i)7o;^ 



7o'iV.(G2-l)Ar2 + A^,72 



Ar2 
9 



1 



(35) 



(36) 



-ai + 2bi 



70 



Equating (1351 ) and (1361 ) the exact MLLA solution of 



1 



reads 

70 



n- 



where 



Si 



Go 



-ai + 2bi 
3 



70 



(37) 



(4ai-/3o). 



Setting 7o/\/^ ^ in (l37l) one recovers the DLA A'^, -independent value ^2^'^ = 4/3. Then, expand- 
ing (|37] ) in the form 1 + 70/ a/ZV^ , one recovers the result from [1 1] for A^s = 1 



where the linear combination of color factors reads 
2 1 2^ 1/55 
'^^ = -r^ + /°+3^^ = 4iVel-9- 



,r^, , 112 ThC^f 32TrC2 



Nr 



+ 



9 A^. A^. 



3 A^. m 



(38) 



(39) 
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3.1.2 Medium-modified Q2 at MLLA and expansion in 0(7o/\/iV7) 

Inserting Tvf ^ = Qs^Vg and using the MLLA expression for the ratio (fTTl ) in (|33] ). it is straightfor 
ward to obtain 

Q2-I 



G2 — 1 Ci? 



(40) 



Expanding (l40l ) in the form 1 + jo/^/N^ and setting A'^g = 1, one recovers the result from [11] in the 
vacuum. Indeed, 



(41) 



where, in agreement with [11], we obtain the combination of color factors 

5 1, 1^ 1/55 ^TrCf STrCI 



(42) 



In Table [3] we display ci together with ci (l39l ) for n/ = 3,4,5. As for the medium-modified MLLA 





Cl 


Cl 


3 


0.485 


0.495 


4 


0.477 


0.491 


5 


0.469 


0.486 



Table 3: Coefficients ci and ci. 



expression r = Ng/Nq ( fTSl ). hard corrections to the MLLA second multiplicity correlators G2 and Q2 
are suppressed by the factor 1 / ^/Wl, while the leading double logarithmic terms (70 / \fN~s —>■ 0) remain 
unchanged and equal the vacuum result 

^2 = 1 + ^, A = g{Cg = N,), A = q{Cg = CF). (43) 

Thus, our MLLA predictions for the medium-modified second multiplicity correlators follow the char- 
acteristics of the jet quenching. 

3.2 Next-to-MLLA evolution equations for the multiplicity correlator 

To obtain the equations we proceed like in paragraph 12.21 and use results from subsection 13.11 Indeed, 
by further pushing the perturbative series, one can improve the account of the energy balance. We 
replace, in the hard splitting region Y ^> Inx ln(l — x), N{Y + Inx) iN{Y + ln(l — x))) and 
iV(2)(y + In 2;) (iV(2)(y ^ in(l - x))) by N(Y) + -^N(Y) Inx... {N{Y) + -^N(Y) ln(l - x) . . .) 
and iV(2)(y) _^ ^]y(2)(Y) Inx... {N^^\Y) + -^M^\Y) ln(l - x) . . .) respectively in the system 
( I30I31I ). while the dependence on Inx is kept on the singular piece Ng/x. After integrating the regular 
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terms over x, the medium-modified NMLLA approximate system of two-coupled evolution equations 
for the gluon and quark multiplicity correlator reads 



+ To ((ai - h) 
d 



ai\^-(i.li)+a',{Ns)^,\NP 



Nf^ - Nf) = ^70^ (at, - ai - Poll) + ^Ns) 



dY^ 



where the term oc ^3 follows from the MLLA result (1401 ). 



Cf 



1 + ( ^ai - 2di 



70 

/n: 



The constants are the following: 



b2{Ns 



2TrCf n 
13 / Cf 

Nr 



2 1 



Cf\ n 



Nr 



(44) 



(45) 



(46) 
(47) 
(48) 



The terms oc ai/3o, a2{Ns), (oi — 6i)/?o, ^2(^5), ^3(^5) in (l44l ) and the ones oc diPo, 02(^5) 
in (|45] ) are 0{jq) corrections which better account for energy conservation. We remind that ~ 
0{{NsasT/^) and that terms oc /?o arise from the running of the coupling constant as{Y). Moreover, 
these constants take A'^s -dependence for the reasons explained in section 



3.2.1 Medium-modified G2 at NMLLA and expansion in O(7o/v'-^s) 

Setting Ng"^^ = G2Ng in (l44l) . the equation can be solved iteratively by making use of (|2TI ). the MLLA 
formula ^ for G2 and the leading DLA limit C^^^ = 4/3; moreover, we expand the series up to terms 
0(0^). The l.h.s. of (|44l) can therefore be written in the form, 



l.h.s. = ^iiC2 - 1) 4iV, - (8c7i + Po)^/Ns^o + 2 (fTi(2ai + /3o) - 4a2(iV,)) 7o 



Ar2 
9 



+ 2/3oCi7o'A^o- 



The r.h.s. reads 
r./i.s. = 



1! 



Ns{G2 -l) + N,-2 Qai + 61^ 7a^7o + Qai + 61^ (2c7i + /?o)7o 



4 , 



Equating (l49l) and (l50l ) we find the new exact NMLLA solution of (|44l) . 



G2 



(49) 



(50) 



(51) 
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where the following combinations of color factors have been written in the form 



(ai - /3o) 0V; 



+ 4Qa'2(iV,) + 62(iV,)) Ns, 
= 2(ai(2ai+/3o)-4a2(iV,)). 



(52) 

(53) 
(54) 



When the MLLA coefficients 6i, 62 and NMLLA 63{Ns) and 6i{Ns) are evaluated in the vacuum (iVg = 
1) for Uf = 3, we find respectively 5i = 2.453, 82 = 2.991, 63(1) = 2.818 and 64(1) = 3.766. In 
particular, 61 ~ ^3(1) while the NMLLA ^4(1) becomes bigger than the MLLA 62- It was shown in 
the KNO problem that MLLA corrections increase like ~ ^^foLg ~ 2) while NMLLA like ~ k^Ug 
(k^ = 4) as the rank of the correlator, which coincides with the number of particles triggered in the 
shower, increases [28]. It may be the reason why sizable NMLLA coefficients are found in this picture. 
Moreover, as the rank k of the correlator increases, 0{^/ag) corrections become of the same order of 
magnitude than the leading DLA and perturbation theory fails. Therefore and in general, MLLA and 
NMLLA corrections for the less inclusive multiplicity correlator of any rank k are more sizable than 
those of the more inclusive average multiplicity. That is the reason for, the exact numerical solution of 
the evolution equations [7,29] becomes interesting. 
Expanding (ISTI ) in 70 / in the form 1 + 70/ + Jq/Ns, "^e obtain 



Go 



1 



ci- 



70 



+ C2{N,)j^ + 



N. 



3/2 



(55) 



where 



C2{Ns 



1 



364{N,) - 36i{Ns)62{Ns) + QS^iNs)) . (56) 

Setting = 1 in (I56ll and taking nj = 3, 4, 5, we recover the values 02(1) = 0.0372, 0.0609, 0.0838 
obtained in the vacuum [13]. Moreover, in (l55l) . the sign of successive terms change as higher order cor- 
rections are added to the series. Consequently, it should be wondered whether this result can drastically 
be affected as higher order terms are incorporated to the series at current energy scales. The highest 
energy scales reached at the LHC and measured by the ALICE and CMS experiments at CERN will 
provide more reliable comparisons with our predictions than current experimental studies at RHIC. 



3.2.2 Medium-modified Q2 at NMLLA and expansion in O (70 / V^s) 

The solution of (1451 ) can also be obtained by setting A'^^-' = Q2Nq in the equation, using (ISTl ) and taking 
the MLLA formula for G2 (l4TI ). one has 

El 

and 

r.h.s. = ^7^ {nsG2 - ^0V^ai7o + 2 0ai(2ai + /3o) + aici + ^a2{Ns)N^ 7^^ N^. 
After equating ( |57] ) and ( |58] ) we obtain the new exact analytical solution of ( |45] ) 



l.h.s. = 2/3o^ci7o'iV^ + 7o(Q2 - 1) (47V. - 82^/^^^ + <54(iVs)7o " 4/3ori7o) 



Q2 



Cf 



Go -6 



70 



^2 ' 



(57) 



(58) 



(59) 
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where (see ( |26b ) 



Moreover, 



— = 1 — ri —= — roily s) — ^ 



5i = —di 
3 



53{Ns) = 2 (^-ai(2ai + /?o) + aici - /3oCi + -a2(iV,)iV, ) , 
.54(A^.) = <54(A^.)-4/3ori, 



(60) 

(61) 
(62) 



and G2 should be taken from (|5T]) . As before, the size of NMLLA coefficients 53{Ns) and 54(A'^5) in the 
vacuum are quite sizable, fornj = 3 one has indeed, 53(1) = 3.598 and ^4(1) = 3.210, which are close 
to k"^ = 4, where k = 2 labels the rank of the second multiplicity correlator. 

Performing the same expansion in 70 / \/iV^ we obtain the result 



Cf 



N.. 



3/2 



(63) 



where the expression for C2{Ns) follows from (l59l ): 

1 



16 



ri 



ci + 2 + -n 
1 



+ -{c2{Ns) + 63{Ns)) + ^(ci + 2) - -{2r2{Ns) - rf 



(64) 
(65) 



Accordingly, setting A'^s = 1 in we find the values in the vacuum £2(1) = 0.215,0.222,0.229 
respectively for nj = 3, 4, 5 like in [13]. The sign of successive terms added to the series ( |63] ) shows the 
wide oscillating property. We give the values of C2(A'^s) and C2{Ns) in Table ID 









C2{Ns) 






C2{Ns) 




3 


-0.258 


- 0.016^ 


/W, + 2.505Ns 


- 2.1937V2 


-0.168 + 0.005^ 


/7vr + 1.962iVs 


- 1.584iV2 


4 


-0.236 


- 0.022\ 


/Ws + 2.513Ns 


- 2.1937V2 


-0.146 + 0.007\ 




- 1.5847V2 


5 


-0.215 


- 0.029\ 


fWs + 2.521A^3 


- 2.193A^2 


-0.126 + 0.009\ 




- 1.584A^2 



Table 4: Coefficients C2{Ns) and C2{Ns). 



3.3 NLO and NNLO results on G2 and Q2 

The MLLA and NMLLA predictions for G2{Q) dSSJl and Q2{Q) dM]) are depicted respectively in Fig.[3] 
and Fig.m At MLLA, the second multiplicity correlator increases as A'^s > 1 and approaches the asymp- 
totic regime A2 = I + Indeed, as for the MLLA ratio r{Ns) (fTSl) . the hard corrections 0(70) are 
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Figure 3: MLLA (1381 ) and NMLLA (l55l) second multiplicity correlator inside a gluon jet as a function of 
Q = EQ in the vacuum {Ng = 1) and in the medium {Ng = 1.6 and Ng = 1.8). 



suppressed by 1/ ^/Wl, such that the production of soft and coUinear hadrons is enhanced, while that 
of hard coUinear hadrons is more restricted. As before, these results provide evidence for the softening 
of jets in the nuclear medium. However, the NMLLA results (I55I63I ) follow the behavior described in 
section [23] for r{Ns). As Ns > 1, the correlators decrease, one finds indeed the rough dependence 
C2{Ns) ~ —Ng, C2{Ns) ~ —N^ (see TablelDl, which in both cases leads to the unavoidable decrease of 
A2 as Ng increases. This result follows from the wide oscillating property of the perturbative series: it is 
wider in the medium than in the vacuum. That is the reason why, the more physical MLLA trends can be 
recovered either by incorporating higher order terms or by numerically solving the evolution equations 
(ISOEB like in [7,29]. 

Another interesting feature of these observables concerns the shape of the curves. They are roughly 
identical and do not prove to depend on the medium parameter Ng. Moreover, there exists evidence for 
a flattening of the slopes as the jet hardness Q = EQ increases for A'^s > 1 (vacuum and medium). 
This kind of scaling behavior is known as the KNO scaling: it was discovered by Polyakov in quantum 
field theory [31] and experimentally confirmed by e+e" measurements [16] for the second and higher 
order multiplicity correlators. This phenomenon implies a jet energy independence of the normalized 
multiphcity correlators, which is not affected by Ng neither at MLLA nor at NMLLA. 



3.4 Role of higher order corrections 

In this paragraph we comment on some progresses that could be carried out beyond the NMLLA approx- 
imation. We take the much simpler example of the gluon to quark average multiplicity ratio and give the 
rough dependence of the third coefficient r^{Ng) that can be added to the series (l26l ) in the form 

with 

n{Ng)'''$^ N^g[a^{Ng)-~a^{Ng)], 

where 

a^{Ng) V -iV.C(3), ~az{Ng) -^^sC(3). (67) 
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Figure 4: MLLA (HTI ) and NMLLA (l63l) second multiplicity correlator inside a quark jet as a function of 
Q = EQ in the vacuum {Ng = 1) and in the medium {Ng = 1.6 and Ng = 1.8). 



These coefficients follow from (I3I4I) by further expanding the perturbative series and therefore, 

-r3TO^crV3(l-^)c(3)>0. (68) 

Thus, replacing — r3(A'^s) ^cx"^ N'^ in (l66l ). the third term changes its sign and therefore, the MLLA 
trends as A'^s > 1 can be recovered. However, the whole calculation requires the implementation of the 
two-loops coupling constant in the solution, and eventually, the inclusion of the time-like sub-leading 
splitting functions in the evolution equations. Nevertheless, as powers of Ng increase for higher or- 
der terms, the perturbative approach fails and the exact numerical solution of the evolution equations 
becomes necessary. 



4 Conclusions 



In this paper we have dealt with the medium-modified average multiplicity and the medium-modified 
second multiplicity correlators in quark and gluon jets. Our calculations are based on the Borghini- 
Wiedemann model [5], which models parton energy loss in a dense nuclear medium. The average 
multiplicity is found, after multiple re-scattering of the relativistic hard parton in the medium, to be 
enhanced by the factor ^/Wl on the exponential leading contribution. The former leads, in particular, 
to the medium-modified anomalous dimension 7med (7 7med ~ \/^7o)- Corrections to the leading 
double logarithmic contribution of the average multiplicity arise from both the MLLA and the NMLLA, 
which better account for the energy balance and for the running of the coupling constant effects as 
in the vacuum. In particular, the NMLLA average multiplicity distribution is softer at NMLLA than at 
MLLA (see Fig.[T]), such that the available phase space for harder coUinear hadronic production becomes 
restricted. The increase of the average multiplicity at NNLO is driven by the factor oc Ng (see (|23])). 

The MLLA scheme provides a more realistic picture of the jet quenching through the study of these 
observables: such is the case of the medium-modified gluon to quark average multiplicity ratio r = 
Ng/Nq. Indeed, hard corrections are suppressed by the extra factor which leads to restriction 

on production of hard partons in quark and gluon jets. Therefore, r approaches its asymptotic DLA 
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limit ro = Nc/Cp = 9/4 when the coherent radiation of soft gluons is enhanced by the medium. The 
amplitude of the oscillating series turns out to be wider in the medium than in the vacuum at all energies. 
Nevertheless, the shapes obtained at MLLA and NMLLA are roughly identical but the series may require 
the incorporation of higher order corrections. Furthermore, in both approaches, the gluon jets are still 
more active than the quark jets in producing secondary particles but these characteristics are related to 
the jet energy dependence of these observables rather than to the sensitivity to the parameter Ns in the 
nuclear medium. 

The second multiplicity correlators in quark and gluon jets in the medium are also computed at MLLA 
and NMLLA. The multiplicity fluctuations of individual events must be larger for quark jets as compared 
to gluon jets just like in the vacuum. The MLLA corrections are suppressed by 1 /\/iVs, such that A2 
approaches the asymptotic DLA regime as A'^g > 1, reproducing the expected physics. In addition, the 
KNO scaling holds at MLLA and NMLLA in heavy-ion collisions, the flattening of the slopes in both 
the vacuum and the medium is roughly reached for the same virtualities Q > 100 GeV of the jet energy. 
As before, the scaling depends on the energy scale Q rather than on the sensitivity to the nuclear factor 
Ns. At NMLLA, the behavior as A'^s > 1 is inverted, but this output can be cured, either by incorporating 
higher order terms to the series or by exactly solving the evolution equations numerically, but this is out 
of the scope of this paper. 

Finally, our results might lead to more accurate prescriptions for the behavior of these observables in 
the presence of the nuclear environment if the treatment of parton energy loss is improved in the future. 
Furthermore, the study of parton energy loss and medium-modified observables would ideally require 
the re-construction of jets in heavy-ion collisions. Of course, the huge background makes this task 
highly delicate. Nevertheless, thanks in particular to important theoretical developments on the jet re- 
constructions algorithms [32] in a high-multiplicity environment, future analysis at the LHC by ALICE 
[33] and CMS [34] look very promising. 



Acknowledgments: I would like to thank B.A. Kniehl for supporting my stay at University of Hamburg, 
as well as S. Albino, F. Arleo and I. Dremin for enlightening discussions and useful comments on the 
manuscript. 



List of Figures 



MLLA ([T3]l and NMLLA (|23]l average multiplicity as a function of Q = EQ in the 
vacuum (Ns = 1) and in the medium (A'^s = 1.6 and A'^s = 1-8) 

MLLA ([TSl l and NMLLA (l26l ) gluon to quark average multiplicity ratio as a function of 

Q = EQ in the vacuum {Ns = 1) and in the medium (Ns = 1.6 and Ns = 1.8) 

MLLA (l38l) and NMLLA (1551 ) second multiplicity correlator inside a gluon jet as a func- 
tion of Q = EQ in the vacuum (Ng = 1) and in the medium (Ns = 1.6 and Ns = 1.8). II6 



18 



4 MLLA (|4T]) and NMLLA ( [63] ) second multiplicity correlator inside a quark jet as a func- 
tion of Q = EQ in the vacuum (Ng = 1) and in the medium (Ng = 1.6 and Ns = 1.8). 



19 



References 

[1] Yu.L. Dokshitzer, V.A. Khoze, A.H. Mueller & S.I. Troyan, Basics of Perturbative QCD, Editions 
Frontieres, Paris (1991). 

[2] K. Adcox et al. (PHENIX CoUab.), Phys. Rev. Lett. 88 (2002) 022301; 
S.S. Adler et al. (PHENIX Collab.), Phys. Rev. Lett. 91 (2003) 072301. 

[3] C. Adler et al. (STAR Collab.), Phys. Rev. Lett 89 (2002) 202301. 

[4] R Arleo, hep-ph/08101193; 

R. Baier, D. Schiff & B. G. Zakharov, Ann. Rev. Nucl. Part. Sci. 50 (2000) 37; 

A. Kovner & U. A. Wiedemann, in Quark Gluon Plasma 3, World Scientific, Singapore, 

hep-ph/0T0415T] 

M. Gyulassy, I. Vitev, X.-N. Wang & B.-W. Zhang. Eucrth/0302077 1 ibid.; 
A. Majumder, J. Phys. G 34 (2007) S377; 

For a more recent review, see also S. Peigne & A.V. Smilga, hep-ph/0810 570p . 
[5] N. Borghini & U.A. Wiedemann, |hep-ph/0506218J 

[6] N. Armesto, L. Cunqueiro, C. Salgado & W.C. Xiang, JHEP 02 (2008) 048. 
[7] S. Sapeta & U.A. Wiedemann, hep-ph/0809.425 1 . 
[8] U.A. Wiedemann, Nucl. Phys. B 588 (2000) 303. 

[9] V. N. Gribov & L. N. Lipatov, Yad. Fiz. 15, 781 (1972) [Sov. J. Nucl. Phys. 15, (1972) 438]; G. 
Altarelh & G. Parisi, Nucl. Phys. B 126, (1977) 298; Yu. L. Dokshitser, Zh. Eksp. Teor. Fiz. 73, 
(1977) 1216; [Sov. Phys. JETP46, (1977) 641]. 

[10] A.H. Mueller, Nucl. Phys. B 241 (1984) 141; Erratum ibid., B 241 (1984) 141. 

[11] E.D. Malaza & B.R. Webber, Phys. Lett. B 149 (1984) 501; E.D. Malaza & B.R. Webber, Nucl. 
Phys. B 267 (1986) 702. 

[12] I.M. Dremin & V.A. Nechitailo, Mod. Phys. Lett. A 9 (1994) 1471; JETP Lett. 58 (1993) 945. 
[13] I.M. Dremin, C.S. Lam & V.A. Nechitailo, Phys. Rev. D 61 (2000) 074020. 
[14] G. Abbiendi et al., [OPAL Collaboration], Phys. Rev. D 69 (2004) 032002. 
[15] J. Abdallah et al. [DELPHI Collaboration], Eur. Phys. J. C 44 (2005) 31 1. 

[16] HRS Coll., Phys. Rev. D 34 (1986) 3304; AMY Coll., Phys. Rev. D 42 (1990) 737; DELPHI Coll., 
Z. Phys. C-Particles and Fields 50 (1991) 185. 

[17] I.M. Dremin & J.W. Gary, Phys. Rep. 349 (2001) 301. 

[18] D. Acosta et al., Phys. Rev. Lett. 94 (2005) 171802. 

[19] A.N. Safonov (for CDF Collaboration), Nucl. Phys. B (Proc. suppl.) 86 (2000) 55. 

20 



[20] R. Perez-Ramos, F. Arleo & B. Machet, Phys. Rev. D 78 (2008) 014019; F. Arleo, R. Perez-Ramos 
& B. Machet, Phys. Rev. Lett. 100 (2008) 052002. 

[21] I.M. Dremin, O.S. Shadrin, J. Phys. G 32 (2006) 963. 

[22] N. Armesto, C. Pajares & P. Quiroga Arias, hep-ph/0809.4428. 

[23] Ya.l. Azimov, Yu.L. Dolcshitzer, V.A. Khoze & S.l. Troian, Z. Phys. C 27 (1985) 65; Yu.L. Dok- 
shitzer, V.A. Khoze & S.I. Troian, J. Phys. G 17 (1991) 1585. 

[24] R. Perez-Ramos, hep-ph/081 1.2418. 

[25] V.A. Khoze & W. Ochs, Int. J. Mod. Phys. A 12 (1997) 2949. 

[26] S. Sapeta & U.A. Wiedemann, Eur. Phys. J. CSS (2008) 293. 

[27] M. Heinz for the STAR Collaboration, nucl-exp/0809.3769. 

[28] Yu.L. Dokshitzer, Phys. Lett. B 30S (1993) 295. 

[29] S. Lupia & W. Ochs, Phys. Lett. B 418 (1998) 214. 

[30] R. Perez-Ramos, JHEP 06 (2006) 019, and references therein; R. Perez-Ramos, JHEP 09 (2006) 
014. 

[31] A.M. Polyakov, Sov. Phys. JETP 32 (1971) 296. 

[32] M. Cacciari & G.R Salam, Phys. Lett. B 641 (2006) 57. 

[33] ALICE collaboration, B. Alessandro et al, J. Phys. G 32 (2006) 1295. 

[34] CMS collaboration, D. d'Enterria (Ed.) et al, J. Phys. G 34 (2007) 2307. 



21 



